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Abstract 
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theorem' states that the past asymptotic behavior of generic type IX solutions is governed by 
Bianchi type I and II vacuum states (Mixmaster attractor). We give a comparatively short 
and self-contained new proof of this theorem. The proof we give is interesting in itself, but 
more importantly it illustrates and emphasizes that type IX is special, and to some extent 
misleading when one considers the broader context of generic models without symmetries. 
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1 Introduction 



Based on work reviewed and developed in [1] and by Rcndall [2], Ringstrom succeeded, some eight, 
nine years ago, to produce the first major theorem about asymptotic dynamics of type IX vacuum 
and orthogonal perfect fluid models [3, 4]: The Bianchi type IX attractor theorem, which states 
that the past asymptotic behavior of generic type IX solutions is governed by Bianchi type I and II 
vacuum states (which constitute the Mixmaster attractor); see Theorem 5.1. In this paper we 
provide a comparatively short and simple new proof of this theorem. Our proof rests on three 
cornerstones: 

(i) We introduce new bounded variables that yield a relatively compact state space, which 
eliminates several of the complications associated with the unbounded variables used by 
Ringstrom. 

(ii) We make systematic use of the Lie contraction hierarchy of invariant subsets admitted by 
Bianchi type IX where monotone functions restrict the asymptotic dynamics to boundaries 
of boundaries. 

(iii) We systematically employ methods, arguments, and results from the theory of dynamical 
systems, see e.g. [1, 5, 6, 7] and references therein. 

However, to find a more succinct argument is not our primary motivation to re-investigate Bianchi 
type IX. Our proof demonstrates that Bianchi type IX is special in comparison with the other 
oscillatory Bianchi models, i.e., type VI_i/g and type VIII. This special nature of Bianchi type IX 
is associated with the geometric condition that the three structure constants have the same sign, 
which results in an extraordinary simplification of the problem (which is easily overlooked) and 
thus makes the treatment of Bianchi type IX models relatively straightforward. The fact that 
Bianchi type IX is rather special has broad ramifications for our understanding of generic spacclikc 
singularities [8, 9]: Bianchi type IX is probably not quite as good a role model for the asymptotic 
behavior of generic inhomogeneous models as is commonly asserted. 

Ringstrom's Bianchi type IX attractor theorem is a remarkable theorem. However, it is imperative 
to point out that this theorem has limited implications, e.g., nothing is rigorously known about 
dynamical chaotic properties (although there are good grounds for beliefs). All claims about chaos 
in Einstein's equations rest on the (plausible) belief that the asymptotic dynamics of Einstein's 
equations is reflected by a discrete map (the Kasner map), but this is far from evident and has 
not been proved so far. These and related issues are discussed in the paper 'Mixmaster: Fact and 
Belief [10]. The present paper, however, concentrates on rigorous results. 

This paper is essentially self-contained. In Section 2 we begin with a brief discussion of the 
dynamical systems approach, where we establish the connection with the metric approach. We 
briefly introduce Hubble-normalized variables, but we concentrate on a set of new bounded variables 
that yield a relatively compact state space — this is the first cornerstone of our analysis. In Section 3 
we discuss the levels of the Bianchi type IX Lie contraction hierarchy: Apart from reviewing results 
on the Bianchi type I and type II subsets, we present a thorough analysis of Bianchi type VIIo- 
The proofs we give are novel and, in particular, independent of results on Bianchi type IX — this 
is the second cornerstone of our analysis. In Section 4 we present the results of the local analysis 
of the fixed points of the dynamical system and discuss non-generic asymptotically self-similar 
behavior. Section 5 is the core of this paper: We present a new and succinct proof of the Bianchi 
type IX attractor theorem (which is stated as Theorem 5.1). The proof is considerably shorter 
than the proof given by Ringstrom [4] (which is in turn based on results in [1, 2, 3]). In the proof 
we systematically employ methods from the theory of dynamical systems — the third cornerstone 
of our analysis. In Section 6 we state and prove a number of consequences of Theorem 5.1. These 
results follow relatively easily from Theorem 5.1 when combined with the knowledge of the flow on 
the Mixmaster attractor (which is the union of the type I and II subsets) . We conclude in Section 7 
with a discussion of the main themes of this paper and we put Bianchi type IX in a broader context; 
in particular we emphasize that the present Bianchi type IX models arc too special in some respects 
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Table 1: The class A Bianchi types are characterized by different signs of the variables (n a , n/},n y ), 
where (a/37) is anv permutation of (123). In addition to the above representations there exist 
equivalent representations associated with an overall change of sign of the variables; e.g., another 
type IX representation is ( ). 



to serve as good role models for generic spacelikc singularities. Throughout this paper we use units 
so that c = 1 and 8irG — 1, where c is the speed of light and G the gravitational constant. 

2 Basic equations 

Consider a vacuum or orthogonal perfect fluid spatially homogeneous (SH) Bianchi type IX model 
(i.e., the fluid 4- velocity is assumed to be orthogonal to the SH symmetry surfaces). It is well 
known that there exists a symmetry-adapted (co-)frame {uj 1 , <1> 2 , uj 3 } such that the metric for 
these models takes the form 

4 g = -dt ®dt + g n {t) uj 1 ® uj 1 + g 22 (t) uj 2 ® uj 2 + g 33 (t) uj 3 ® uj 3 , (la) 
where duj 1 — — h\ uj 2 A uj 3 , duj 2 = — h 2 uj 3 A uj 1 , duj 3 — —h 3 li'Au 2 , (lb) 

and where h a — +1 Va; see e.g. [1, 4] and references therein. 

Remark. Class A Bianchi models of different Bianchi types are characterized by different structure 
constants n a ; the different cases are listed in Table 1. 

Let 

m(t) := m _— , n 2 {t) := n 2 — == , n 3 (t) := n 3 -== , (2) 
ydet g vdet g ydet g 

where dct g = gng 22 g 33 . Furthermore, define 

= -trfc and a a = -k a p + | tr k 8 a = diag((Ti, a 2 , cr 3 ) (=> a °> = °) ' ( 3 ) 

where k a p denotes the second fundamental form associated with (1) of the SH hypersurfaces 
t = const. The quantities 6 and a a p can be interpreted as the expansion and the shear, respectively, 
of the normal congruence of the SH hypersurfaces. The spatial volume density changes according 
to according to dy/dctg/dt = 6*Vdet g. In a cosmological context it is customary to replace 
by the Hubble variable H = 6/3 = — trfc/3. Evidently, for the present Bianchi type IX models 
there is a one-to-one correspondence between the 'orthonormal frame variables' (H,a a ,n a ) (with 
J2 a a a =0) and (g a p, k a(3 ). 

Remark. The orthonormal frame variables (H, cr a ,n a ) can be used to describe any model of the 
family of Bianchi class A models; see Table 1. In Bianchi type VIII and IX, the metric is obtained 
from (ni,n 2 ,n 3 ) by (2); for the lower Bianchi types, the other frame variables, i.e., {H,<j a ), are 
needed as well to reconstruct the metric; see [11] for a group theoretical approach. 

In the perfect fluid case we assume an orthogonal perfect fluid with density p and pressure p that 
satisfy a linear equation of state p — wp. We require the energy conditions (weak/strong/dominant) 



3 



to hold, hence p > and 1 

-l<w<l; (4) 

we exclude the special cases w = — ^ and w = 1, where the energy conditions are only marginally 
satisfied. 2 



2.1 Hubble-normalized variables and equations 

In the Hubble-normalized dynamical systems approach we define dimensionless orthonormal frame 
variables according to 

{■E a ,N a ) = (a a ,n a )/H , n = p/{'3H 2 ). (5) 

Remark. For all class A models except type IX the Gauss constraint guarantees that H remains 
positive if it is positive initially. In Bianchi type IX, however, it is known from a theorem by Lin 
and Wald [14] that all type IX vacuum and orthogonal perfect fluid models with w > first expand 
(H > 0), reach a point of maximum expansion (H = 0), and then recollapse (H < 0). 3 Therefore, 
although the variable transformation (5) breaks down at the point of maximum expansion, the 
variables (£ Q , A Q ) correctly describe the dynamics in the expanding phase. 

Since the past singularity is of particular interest in our considerations, it is convenient to choose 
the time direction towards the past singularity. We introduce a new dimensionless time variable 
t_ according to 

dr^/dt = -H . (6) 

The Einstein fields equations can be reformulated in terms of the Hubble variable H and the 
Hubble-normalized variables [1]. Since H is the only variable that carries dimension, the equation 
for H, 

dH/dr- = (1 + q)H , (7) 
decouples from the equations for £ Q and N a , which are given by the system (as follows from [1]) 

dE a /dr- = (2 -q)E a + 3 S a , (8a) 

dN a /dT- = — (q + 2E Q ) N a (no sum over a) , (8b) 

where 

q = 2E 2 + i(l + 3w)n , E 2 = |(S 2 + S 2 + S 2 ), (9a) 

and 3 S a = i [N a (2N a - N p - Ag - (N - iV 7 ) 2 ] , (a/87) G {(123), (231), (312)} . (9b) 

It is straightforward to show that p oc exp (3[1 + u>]t„) and t_ — > 00 towards the past singularity, 
so that p — > 00 in this limit; see [10]. 

Apart from the trivial constraint Si + £2 + S3 = 0, 4 there exists the Gauss constraint 



Nf + 7V| + iVf - 2(^1^2 + AWs + N 3 N-_ 



l-fi<l. (10) 



Since this constraint can be used to globally solve for fi, the reduced state space is given as the 
space of all (Ei,E2,£ 3 ) and (N 1 ,N 2 ,N 3 ) satisfying (10) and Y, 1 + E 2 + E 3 = 0. The Gauss 
constraint (10) reveals a serious disadvantage of the Hubble-normalized variables: The range of 
the variables is unbounded, i.e., the state space is not relatively compact. 



1 Note, however, that the well-posedness of the Einstein equations (for solutions without symmetry) has been 
questioned in the case — l/3<to<0, see [12]. 

2 The case to = 1 is known as the stiff fluid case, for which the speed of sound is equal to the speed of light. The 
asymptotic dynamics of stiff fluid solutions is simpler than the oscillatory behavior characterizing the models with 
range — i < w < 1, and well understood [4, 13]. (In the terminology introduced below, the stiff fluid models are 
asymptotically self-similar.) We will therefore refrain from discussing the stiff fluid case in this paper. 

3 In the locally rotationally symmetric case it has been proved that the range of w can be extended to 10 > - 
see [15]. There are good reasons to believe that the assumption of local rotational symmetry is superfluous, but this 
has not been established yet. 

4 It is common to globally solve Ei + £2 + S3 = by introducing new variables according to Si = — 2E+, 
S2 = £~f — v3S_, S3 = + \/3X!_, which yields S 2 = E 2 ^ + However, since this breaks the permutation 
symmetry of the three spatial axes (exhibited by type IX models), we choose to retain the variables Si, £2, S3. 
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2.2 Bounded variables 



The preparatory step in our approach to Ringstrom's Bianchi type IX attractor theorem is to 
reformulate the Einstein equations for class A Bianchi models in terms of variables that span 
a relatively compact state space. (The aim is to avoid the problems that are caused by the 
unboundedness of the Hubble- normalized variables in the dynamical system (8); note, however, 
that the key arguments in our proof, e.g., Lemma 5.3, are independent of the choice of variables.) 
Define 

D := \Jh 2 + |(nin 2 + nin 3 + n 2 n 3 ) . (11) 

The Gauss constraint reads ZD 2 = p+ |(crf + a 2 + erf ) + \{n\ + n 2 + n|), hence D > 0. This makes 
it possible to introduce variables that are normalized w.r.t. the 'dominant' variable D instead of 
H, i.e., 

(S a ,N at H) = {tT a ,n a ,H)/D, fi = p/{3D 2 ) . (12) 

By construction, Si + £2 + £3 = 0. 

Remark. Contrary to the Hubble-normalized variables, the new variables (12) are globally well- 
defined. 

We choose a dimensionless time variable f_ , which is defined by 
accordingly, f_ is directed towards the past singularity. 

The quantity D decouples from the other equations for dimensional reasons, 

J^- = [(1 + q)H + F]D. (14) 

The remaining dimensionless system of equations reads 
dH 

-—- = q(l - H 2 ) — FH , (15a) 

GtT_ 

S a [(2 - q)H -F]+ 3 S a , (15b) 



dN a 



-N a [qH + 2E Q + F] , (no sum over a) (15c) 



where 



q = 2t 2 + i(l + 3w)n , t 2 = \{t\ + t 2 + t 2 ) , (16a) 

F = i (A>iA> 2 S 3 + N x t 2 N 3 + SiA^iVs) , (16b) 
3 S a = i [N a (2N a -N - A> 7 ) - (N p - A> 7 ) 2 ] , (api) e {(123), (231), (312)} . (16c) 

The vacuum case is characterized by f2 = while f2 > in the fluid case. Note that the 'deceleration 
parameter' q is non-negative. 

In contrast to the system for the Hubble-normalized variables, there exist two non-trivial con- 
straints, the Gauss constraint and a constraint resulting directly from (11), 

^ 2 + T2 (^1 + ^2 +Nl)+ti = l, (17a) 



H 2 + i (NxN 2 + N1N3 + N 2 N 3 ) = 1 . (17b) 

Evidently, in Bianchi type IX, the range of the new variables (12) is bounded. This is true for the 
entire class A: We minimize the expression NiN 2 + NiN 3 + N 2 N 3 appearing in (17b) under the 
condition + iVf + iVf < 12 resulting from (17a); this leads to + + N 2 N 3 > -6 and 
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hence H 2 < 2 for all class A models; restriction to Bianchi type IX yields H 2 < 1, since N a > Va. 
Therefore, the system (15) is a system on a relatively compact state space for the entire class A. 

The dimcnsionlcss state space of the Bianchi type IX orthogonal perfect fluid models with a linear 
equation of state is 5-dimcnsional while the state space of the vacuum models is 4-dimensional. 
The same is true for Bianchi type VIII, while the state spaces of the remaining class A Bianchi 
models have less degrees of freedom; see Table 2. Once the dynamics in the dimcnsionless state 
space is understood, D is obtained from a quadrature by integrating (14), and subsequently the 
metric can be reconstructed. Note that we choose to use (17a) to solve for f2 when f2 > 0; however, 
it is still useful to consider the evolution equation 

dfl/df- = -Q[2qH- (l + 3w)H + 2F] . (18) 

Remark. Since q is bounded as r_ — + oo, the asymptotics of H can be bounded by exponential 
functions from above and below. It follows that the equation dt/df_ = —D^ 1 can be integrated 
to yield t as a function of f_ such that t — > as f _ — > +oo. 

In all Bianchi types except type IX the constraints force H to have a sign for all f_, e.g., H > 
(which entails that these models are forever expanding 5 ). In Bianchi type IX the subset H > 
is a future invariant subset of the state space, which follows from the inequality dH/df^ > on 
H = 0. (There is a close relationship between this fact and the results of [14].) 

Let (a/37) denote any permutation of (123). We define the Bianchi type IX state space B\x as the 
future-invariant set 

B IX = {N a > 0, Np > 0, A> 7 > 0, H > 0, E„, T,p, E 7 , O > 0} , (19) 

where the variables are subject to the constraints (17) and E1 + E2 + E3 = 0; the set £?ix comprises 
both the vacuum subset and the fluid subset, i.e., f2 > 0. 

Setting one or more of these variables (N a , Np, iV 7 ) to zero (which corresponds to Lie contrac- 
tions [16]) yields invariant boundary subsets which represent more special Bianchi types. In this 
spirit, the boundary 6 dB\x of the Bianchi type IX state space is given by 

dB IX = H„U Byiio , (20) 

where 

H = {N a >0,N t3 >0,N 1 >0,H = 0,n>0}, (21a) 

Byno = {N a >O,N >O,Ny = O,S>O,n>O}, (a/3j) e {(123), (231), (312)} . (21b) 

Note that Byu denotes the collection of the three equivalent Bianchi type VII subspaces; if we 
want to refer to one of the subsets in particular we use the notation S^a^ f° r the set {N a > 
0, N/3 > 0, N 7 = 0, H > 0}. The notation is such that the subscript denotes the non-zero variables 
among {N a , Np,N^}; accordingly, B V ii = u ®n 2 n 3 u % 3 JVi- 

The constraints (17) imply that the boundary of the set Ho consists of the three points 

Z Q : H = 0, S Q = 0, N a = 0, Np = A> 7 = VE, Cl = 0, (a/3 7 ) G {(123), (231), (312)} ; (22) 

the boundary of Bvn is the union of the closures of the Bianchi type II subspaces, which are 
given by Bjj a = {N a > 0, Np = 0, iV 7 = 0, > 0}, a = 1, 2, 3, collectively denoted by B u , and the 
three points Z Q , i.e., 

&Bvii =5iiU{Zi,Z 2 ,Z 3 }. (23) 

For completeness, we note that B\ is the set {N a = 0, Np = 0, iV 7 = 0, O > 0}; we have dB\\ = Bj. 
A Bianchi subset contraction diagram for type IX is given in Figure 1. 

5 This excludes the Bianchi type I and type VIIo representations of Minkowski spacctime. 

6 In our context, dB t denotes B*\B t \ the word 'boundary' is chosen in lack of better terminology. Accordingly, 
since B IX = {N a > 0, Np > 0, 7V 7 > 0, H > 0, fi > 0}, 8B 1X is the set where one of the variables (JVq, , Np , A? 7 ) or 
H is set to zero, while the vacuum subset {N a > 0, Np > 0, Nj > 0, H > 0, Q = 0} does not appear in 9Bjx- This 
convention is adapted to the formulation of Lemma 5.2. 
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Type Symbol Range of (N a , Np , iV 7 ) 



I 


Bi 


N a 


II 


B u 


N a 


VIo 


Bvi 


N a 


VIIo 


-Bviio 


N a 


VIII 


Bym 


N a 


IX 


Bix 


Na 



0, Np = 0, N y = 
0, N p = 0, iV 7 > 
0, Np < 0, 7V 7 > 
0, iV/3 > 0, iV 7 > 



State space properties 



H = 1, E 2 < 1 
= 1, E 2 -)-^ 2 < 1 
H € (0, 2), E 2 + i[7V^ - iV 7 ] 2 < i? 2 
5 6(0,1), S 2 + ^[^ 

# G (0,2), E 2 < #_ 2 

e (-1,1), e 2 ^ if 2 



iV 7 ] 2 < H 2 



D 



Table 2: The dimcnsionlcss state spaces associated with class A Bianchi models; here, (a/37) i s 
any permutation of (123). In addition to the above representations there exist equivalent represen- 
tations associated with an overall change of sign of the variables (Ni, N2, N3) or H. The quantity 
D denotes the dimension of the state space (in the fluid case); the dimensionality of the state 
space in the vacuum cases is given by D — 1. The Bianchi type IX state space decomposes into 
a future- invariant half and a past-invariant half. By Bix we denote the future-invariant half, i.e., 
H > 0; cf. (19). 



Remark. The constraints imply that H = 1 on B\ U Bu; hence the dominant variables reduce to 
the standard Hubble-normalized variables for type I and II. 

Remark. Each of the sets B\, Bu, Byu , Bix (and Hq) decomposes into a vacuum subset (f2 = 0) 
and a fluid subset (O > 0). We refer to these subsets using the superscripts vac - and fl \ 

Remark. Note that there does not appear a subset Byi as a boundary subset of .Bix- In the 
concluding remarks, Section 7, we argue that this simple fact has far-reaching consequences. 



3 The boundaries of the state space: Bi, Bu, -Bvn 

In the analysis of a dynamical system that is defined on a (relatively) compact state space, the 
analysis of the boundary subsets is crucial. Often there exists a hierarchy of boundary subsets that 
governs the asymptotic behavior of solutions, see e.g. [17, 18] for examples how such hierarchies 
can be exploited. In the case of Bianchi type IX models, the 'Lie contraction hierarchy' of the 
boundaries is depicted in Figure 1. 



3.1 The Bianchi type I subset 

The Bianchi type I subset Bi is given by N x = 0, N 2 = 0, N 3 = 0, hence H = 1 and f2 = 1-S 2 > 0. 
The vacuum subset, O = 0, consists of a circle of fixed points — the Kasner circle K°, which is 
characterized by S 2 = 1. Each fixed point on K° represents a Kasner solution (Kasner metric). 
There exist six special points: Q Q are given by (S a ,,S i g,E 7 ) = (—2,1,1); the Taub points T a 
are given by (E a ,E^,E 7 ) = (2,-1,-1). The former are associated with locally rotationally 
symmetric (LRS) solutions whose intrinsic geometry is non-flat; the latter correspond to the flat 
LRS solutions — the Taub representation of Minkowski spacctimc. 

The Bianchi type I perfect fluid subset is the set 1 — fl = E 2 < 1. From (18) it is straightforward 
to deduce that there exists a central fixed point, the Friedmann fixed point F, given by Ei = 
E2 = E3 = 0, which corresponds to the isotropic Friedmann-Robertson- Walker (FRW) solution. 
Solutions with < E 2 < 1 are given by radial straight lines originating from F and ending at K°. 
These results rely on the assumption w < 1, see (4). 
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D A 




Figure 1: Subset contraction diagram for Bianchi type IX. D denotes the dimension of the dimen- 
sionless state space for the various models with an orthogonal perfect fluid with linear equation 
of state; the associated vacuum subsets have one dimension less; see also Table 2. The notation 
is such that the subscript of B* denotes the non-zero variables, e.g. SjviJV 2 denotes the type VIIo 
subset with Ni > 0, N 2 > and N 3 = 0. 




Figure 2: Projection of the flow on the type II subsets B™-, B v ^\ B v ^ c - onto (Si, S 2 , S 3 )-space. 



3.2 The Bianchi type II subset 

The Bianchi type II subset B\\ has three equivalent representations: B^ !l , B^ , B^ 3 . Let us 
consider , which is given by N a = Np = 0, iV 7 > (hence H = 1). On Bft , the Gauss 
constraint E 2 + yj-/V 2 +0 = 1 can be used to replace N 7 by as a dependent variable. The 
system (8) thus becomes 

^. = (2-q)t a/p + a S a/p , ^ = (2-q)^+%, ^- = -n[2q-(l + 3w)}, (24) 

where q = 2£ 2 + ±(l + 3w)0 and 3 S a/p = -4(l-E 2 -0), 3 S* 7 = 8(l-E 2 -0); we have £ 2 + < 1. 

Let us first consider the vacuum subset B™ c ', i.e., = 0. There do not exist any fixed points on 
£?™ c ', but the boundary coincides with the Kasner circle K°. The orbits of (24) form a family 
of straight lines in 23™ c '; each orbit is a heteroclinic orbit, since it connects two different fixed 
points. If the initial point is Q 7 , the final point is T 7 (LRS orbit); the points T a and Tg are not 
connected with any other fixed point (they are 'fixed points' under the present 'vacuum Bianchi 
type II map'), see Figure 2. 

While there do not exist any fixed points in $™ c ', there exists one fixed point in B^ with O > 0, 
the Collins-Stewart fixed point CS 7 , which corresponds to one representation of the LRS solutions 
found by Collins and Stewart [19]. CS 7 is given by (E a , Eg, E 7 ) = i(l + 3u>)(l, 1, —2) and = 
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1 — + 3w) (which yields iV 7 = |V1 — w\/l + 3w ). The fixed point CS 7 is the source for all 
orbits in . For a detailed discussion of these results see [1]. 

3.3 Bianchi type VII : A new analysis 

A detailed analysis of the dynamics of Bianchi type VIIo models is essential for an understanding 
of type IX asymptotic dynamics. To underline the importance of Bianchi type VIIo, we present a 
new analysis of the global dynamics of type VIIo solutions; the proofs we give are novel and, in 
particular, independent of results on Bianchi type IX. (The proof given in [4, 20] relies on results 
on Bianchi type IX.) 

The dynamical system on -Bvn (— ^JViJV ^ ^n 2 n U N t ) does not admit fixed points with 
Cl > (while, of course, there exist fixed points on <9_Bvn )- However, each of the three vacuum 
subsets contains a line of fixed points with N a = Np > 0: 

TL 7 = jo < H < 1, E a = T,p = -H, £ 7 = 2H, N a = Np = vWl - H 2 : N 1 = 0, fl = o| . (25a) 

These lines of fixed points ('Taub lines') connect the Taub points T 7 on the Kasner circle with the 
fixed points Z 7 . Each of the fixed points on TL 7 represents the Minkowski spacetime in a Bianchi 
type VIIo LRS symmetry foliation; this is in analogy to the Taub points T 7 themselves. (Note that 
T, a = Tp = —H, £ 7 = 2H corresponds to S Q = £3 = —1, S 7 = 2 in Hubble-normalized variables. 
N a = Np G (0, v6) corresponds to N a = Np G (0, oo).) The second family of LRS vacuum subsets 
are the three sets 

QC-y = jo < H < 1, S a = tip = H, S 7 = -2H, N a = Np = Vq^I-H 2 , N 7 = 0, Cl = o} , (25b) 

which connect the exceptional points Z 7 with the non-flat LRS fixed points Q 7 . These lines 
are not lines of fixed points, since dH /df_ = AH 2 (1 — H 2 ), but heteroclinic orbits Z 7 — > Q 7 . 
The corresponding solutions represent the non-flat LRS Kasner solutions in a Bianchi type VIIo 
symmetry foliation. 

The vacuum case f2 = 0. On the representative S™ c ^ of -E?vii cons ider the function (2H — 

S 7 ). This function is non-negative, since S 7 G [—2H,2H], which follows from the constraints 
when we use the relation 12S 2 = 3S 2 + (S Q — Tp) 2 ; furthermore, (2H — S 7 ) is increasing on 

F£ C ^\(TL 7 UK°), since 

J- (2H - S 7 ) = A. (2H - S 7 ) [(N a + Npf (2H - £ 7 ) + (N a - N ) 2 (2H + S 7 )] . (26) 

The monotonicity principle 7 implies that the a-limit is a fixed point on TL 7 = {T 7 } UTL 7 U {Z 7 }, 
since the fixed points on the Kasner circle are excluded as possible a-limit points because they are 
(transversally hyperbolic) saddles and sinks. 

Remark. In addition, the monotonicity principle implies that the w-limit of every orbit in Bff C ff 

(except for the fixed points on TL 7 themselves) is a fixed point on the Kasner circle (where T 7 is 
excluded) ; clearly, only the transversally hyperbolic sinks come into question. 

The results concerning the a-limits of orbits that arc obtained from the monotonicity principle can 
be strengthened considerably: 

Lemma 3.1. The only type VIIo vacuum orbit that converges to Z 7 as f_ — > — oo is the orbit 
QC 7 ; any other orbit converges to one of the points on TL 7 as f_ — > —oo; conversely, each point 
on TL 7 is the a-limit set for a one-parameter set of orbits. 

7 The monotonicity principle [1] gives information about the global asymptotic behavior of solutions of a dynamical 
system. If M is a C 1 function on the state space X that is strictly decreasing along orbits, then 

w(as) C {£ £ X\X | lim Af(C) + supM} , a(x) C {§ £ X\X \ lim M(Q + inf M} 

C—>€ x C— »€ X 

for all x £ X. 
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The proof of these statements requires a more detailed analysis of the dynamical system (in the 
neighborhood of TL 7 , in particular) . The reader who is not interested in these details might prefer 
to continue with the discussion of the Bianchi type VIIo fluid case below; note, however, that the 
proof we give here is independent of the proof given in [4, 20] (and, in particular, it is completely 
independent of results on Bianchi type IX). 

Proof. Let us introduce a set of variables that are adapted to the special features of the Bianchi 
type VIIo state space, using the Hubble-normalized formulation as a starting point. 8 Consider 
B}? C m and let 

S Q - Y, p = 2\/3 sin <d cos 2tp, £ 7 = 2 cos «?, N a - Np = Vl2 sin # sin 2ip, N a Np = 3( 2 , (27) 

where (■&, ip, £) £ [0, n] x [0, n) x (0, oo). The condition i9 = tt yields Q£ 7 , while •& = yields TL 7 . 
The constraint is automatically satisfied by this choice of variables; the dynamical system takes 
the form 

d-d 

- — = (l-cosV0(l -cosi?)sintf, (28a) 



d ^ „ ■ a L( C 2 , 1-COS^ 

4smiM3 



— ~X Dill if \l *J\ n 

ar_ V Vsin $ 

dC 
dr. 



(C 1 — COS ID \ 
— V- + — -) -2sin^(l-cosi9) , (28b) 
sin v 2 / 



= -C [2(1 - cos t?) - sin 2 i? (1 - cos V)] ■ (28c) 



Consider a non-LRS orbit in 23™°^, i-e., assume < d < tt. The function £/sini9 appearing 
in (28b) is a monotone function, 

d C C 

— ^i-=-^i-(l + coB^(l-cosi?); (29) 
dr_ sm» suit? 

in particular, for sufficiently small t_, £/sin?9 is bounded from below by a positive constant. The 
function -d is monotonically decreasing in the reversed direction of time; most importantly, $ — ► 
as t_ — > — oo. (Proof: Assume the contrary, i.e., d — > = const > as r_ — > — oo. 9 The 
expression 

d r. C ^ N 



— log log tan- =-2(1 -cost?) (30) 

g(t_ V sin # 2/ 

converges to the limit c = —2(1 — cos $00) as t_ — > —00. We thus obtain the asymptotic behavior 
( ~ e _CT ~ as r_ — > —00; accordingly, from (28b), dip/dr- ~ — e~ cr ~ and hence -0 ~ e _CT ~ 
as r_ — » —00. Therefore, the integral of cosip behaves asymptotically like the cosine integral 
Ci(e -CT ~); in particular, the limit exists as t_ — ► —00. This, however, contradicts the assumption 
$ — * $00 > because of (28a).) Likewise, for sufficiently small r_, the function ip is strictly 
monotonically decreasing, i.e., dip/dr^ < 0. This is a consequence of the property z? — > as 
t_ — > —00 (and C/sini? being bounded away from zero). 



In the limit — ^ 0, the system (28) takes the form 
-^ = i(l-cos^ 3 (l + 0(tf 2 )), 



— - ±(l-cos^)tf 3 (l + 0(tf 2 )) , A- ^-([cos^^+O^ 4 )] , (31a) 



rfr_ V Vsin 2 ?? 



W3(-^ + ^) + ( ^). (31b) 



We introduce an alternative time variable, c_, by requiring der_ = i9 2 gJt_. Evidently, er_ is a 
monotone function of t_; most importantly, cr_ — > —00 as t_ — > —00. (Proof: Assume the 



8 It follows from Section 2 that = T,g/H and JV5 = Ng/H, 8 = 1,2,3; in particular, S 7 
time variables satisfy dr_ = Hdf— . 

9 This assumption does not suffice to conclude that dD/dr— — > as t_ — > — 00, since the second derivatives are 
in general not bounded, nevertheless, the heuristic reasoning is correct: d-S/dr^ approaches zero as t_ — » —00 and 
thus t\> approaches a multiple of 27r in this limit, cf. (28a), which is a contradiction to (28b). 
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contrary, i.e., <7_ — > (Too > — oo as t_ — > — oo. Expressed in tr_, Eq. (31a) reads d log $/der_ = 
(1 — cos')/')(l + 0(i9 2 )). However, this contradicts the fact that $ — > as ct_ — > (Too.) Using the 
variable <r_ (which can be identified as a logarithmic time variable), the system (31) becomes 

dh£± = !(l-co S V0(l + O(tf 2 )), li^ = _ C o S V; + 0(,9 2 ), (32a) 

CvC/ tXC/ — 



der_ V V s in 2 

in addition, 

^( l0g ^)=4( 1 + C ° S ^(l + ^ 2 )). (33) 

The function C/sin$ is monotone, cf. (29); furthermore, C/sin?9 — > oo as er_ — > — oo. (Proof: 
Assume the contrary, i.e., C/sint9 — > eioo = const > as cr_ — > — oo. Expressing (30) in terms of 
er_ and integrating this equation yields ~ e CT ~ as er_ — > — oo. Consistency with the differential 
equation for cf. (32a), requires the integral of cos'0 as <r_ — > — oo to converge to — 1. This, 
however, contradicts the differential equation for ip, cf. (32b).) 



' =-^-^(1 + 0(1)), (34) 



Finally, we obtain 

gJ<7_ sin $ 

where the function on the r.h.s. goes to — oo monotonically as c_ — » — oo (since and C/sint? 
are monotone and grow beyond all bounds as cr_ — ► — oo). We conclude that the integral 



dcr'_ cosip(a'_) (35) 

exists. (A statement like this can be regarded as a continuous version of the Leibniz criterion 
for alternating series.) Therefore, by integrating the differential equation for C it follows that C 
converges to a positive constant as cr_ — ► — oo. (The solution of the equation for i? implies that 
the lower order terms cannot contribute.) 

We therefore obtain that the Hubble-normalized variables N a and Np converge to (one and the 
same) constant as r_ — > — oo; consequently, the a-limit set of a non-LRS orbit is a point with 
E a = = — 1, E 7 = 2, N a = N/3 > (and iV 7 = 0). Expressed in terms of the dominant 
variables, the a-limit set of a non-LRS orbit in BJ^ft is a point on TL 7 . □ 

Remark. Based on this information it is straightforward to consider and analyze the asymptotic 
system of differential equations, which arises by inserting the asymptotic behavior of solutions 
into (28). The asymptotic oscillations of solutions can be read off directly. 

The fluid case O > 0. On the representative Bft ft of -Byiio we adapt a monotone function 
found by Uggla, given in [1], 

Co = (2H - vt, 1 )~ 2{1+v] n (N a N p ) v with v = \(l + 3w). (36) 

Note that < v < 1 since -| < w < 1; hence 2H - u£ 7 > if H > (since E 7 € (-2H,2H) 
because of the constraints) and thus < Co < oo on B ft a ft \{Z-y} ■ On this subset we find 



df 3 



^-Co = - (2H - vS 7 ) 'Co 
= - (2H - vt^y 1 Co 

d( /df- =0 



|(1 - v 2 )(t a - t p ) 2 + 2 (2vH - E 7 



|(1 - , 2 )(A 2 - TV 2 ) 2 + (1 - vf (8«fi - l(N a - NpfY 



hence Co is monotonically decreasing on Bft ^ . 

This allows us to apply the monotonicity principle: Using that Co — on B\\ and on the vacuum 
subset B^ c jy^ we conclude that the a-limit of every orbit in Bft ft is one of the fixed points Z 7 . 
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The asymptotic approach to the fixed points Z 7 provides an example for asymptotic self-similarity 
breaking, see [21] for details and explicit decay rates. (Note that the points Z 7 do not appear as 
fixed points in the Hubble-normalized approach; instead, they are associated with 'infinity' in the 
Hubble-normalized state space and thus do not correspond to self-similar solutions). 

The results on the w-limit sets of type VIIo orbits are not needed in the proof of the Bianchi type IX 
attractor theorem below. However, for completeness we now analyze the possible w-limit sets of 
orbits in Bj^ ^ . The monotonicity principle, applied to the function (36), implies that the cj-limit 

sets of orbits must be contained in Bfj a or B^ fj (which are part of £?n) or/and in the vacuum 
subset . However, on only TL Q is admissible. To prove this assume that there exists 

an orbit 7 that has an w-Umit point P on B^°j^ and that this point does not lie on TL 7 . Then the 
orbit through P and the w-limit set of that orbit must be contained in u>("j). As proved above, this 
w-limit set consists of one single fixed point Kp on the Kasner circle (which acts as a sink, when 
viewed as a fixed point on the [closure of the] vacuum subset). Since Qr 1 dO./df- = —3(1 — w) 
on the Kasner circle, Kp is a (transversally hyperbolic) sink also when viewed as a fixed point on 
Bfj a j] g . The conclusion that the sink Kp is in ^(7) contradicts the assumption P g ^(7). This 
leaves only the fixed points on TL 7 as possible w-limits on the vacuum subset B^°^ . In fact, for 

each fixed point L 7 £ TL 7 there exists exactly one orbit in B^ ^ converging to L 7 as f_ — » 00. 
This follows from the center manifold theorem by noting that the fixed point L 7 acts as a center 
saddle (which is because ti~ 1 dti/d,T- = —3(1 — w)H along TL 7 ). One easily identifies the orbit 
that converges to L 7 as being an LRS orbit (the LRS case is exactly solvable). We have thus shown 
that convergence to TL 7 occurs for a non-generic set of orbits; the w-limit set of every other orbit 
lies on the sets B^ a and Bft (and is thus of type I or II). 

The details are as follows: The remaining LRS solutions either converge to Q 7 (a one-parameter 
family) or converge to F as f_ — > 00 (one solution). The local analysis of these fixed points, 
see Section 4 for the general case, implies that the former are embedded into a two-parameter 
families of orbits converging to Q 7 , while the latter are embedded into a one-parameter families of 
orbits converging to F. Furthermore, there exist two (equivalent) orbits converging to each of the 
fixed points CS Q , CS/3. 10 The generic scenario, however, is convergence to one of the transversally 
hyperbolic sinks on the Kasner circle. To prove this statement suppose that there exists an orbit 7 
that possesses an w-limit point P on B^ a (B^ ) with P 7^ CSq, (P 7^ CS/3). Then the orbit through 
P and the a-limit set of that orbit must be contained in w("f). As stated in Section 3.2, this a-limit 
set coincides with CS a (CS/3) in &N a (^n^)- As this point is a hyperbolic saddle in Bfj a fq i}1 the 
orbit converging to CS Q (CS/3) as f _ — > 00 must be contained in 0^(7) as well. However, this is 
a contradiction to the fact that ^(7) is a subset of Bff a U B^^. Analogously, one can prove that 
^(7) and the interior of B\ ■ are disjoint, which leads to the statement. Summing up, the cj-limit 
of a generic orbit in -Byii * s onc or " ^ ne transversally hyperbolic sinks on the Kasner circle. 



4 Non-generic solutions: Asymptotic self-similarity 

In the previous section we have identified the fixed points associated with the system (15) on the 
Bianchi boundary subsets of Six- A local dynamical systems analysis of the fixed points shows 
whether or not these points attract type IX orbits in the limit f_ — > 00. 

K° Evaluated on the Kasner circle, Eq. (15) implies N~ l dN a /df- |j<o = —2(1 + S Q ) (a = 
1,2,3) and fi _1 dfi/df_ | K o = -3(1 - 10). Each fixed point K on K°\{Ti, T 2 , T 3 } is a 
transversally hyperbolic saddle that has one unstable mode and three stable modes. The 
unstable manifold of K coincides with a vacuum type II orbit, see Figure 2; the three- 
dimensional stable manifold is contained in Byu a . (The one-dimensional center manifold 
is K° itself.) Therefore, there do not exist any type IX solutions that converge to K as 
f_ — > 00. The Taub points {Ti,T2,T3} are not transversally hyperbolic. Each Taub point 

10 Equivalence of orbits refers to the discrete symmetries of the problem that are associated with interchanging 
the axes. 
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T a possesses a two-dimensional stable manifold and a three-dimensional center manifold. The 
(closure of the) two-dimensional stable manifold of T a coincides with the LRS subset of Bft a 
(which contains the two special orbits Q Q — ► T Q and F — > T Q on the vacuum subset of Bft 
and on the Bianchi type I fluid subset, respectively.) The three-dimensional center manifold 
coincides with the vacuum subset of Bft ft (which is a vacuum Byii set). Therefore, the 
center manifold reduction theorem [6] reduces the problem to analyzing Bianchi type VIIo 
vacuum dynamics. In Section 3.3 we have shown that T Q is excluded as an w-limit point 
for orbits in ; the existence of the monotone function (26) implies that T a is a center 

saddle in Bix- Consequently, there do not exist any type IX solutions that converge to any 
of the Taub points as f_ — > oo. 

F Eq. (8) implies that E^dEa/df- | F = |(1 - w) and N^dNa/df- | F = -|(1 + 3w). There- 
fore, F is a hyperbolic saddle that possesses a two-dimensional unstable manifold, which 
coincides with the Bianchi type I subset, and a three-dimensional stable manifold. Accord- 
ingly, F attracts a two-parametric family of type IX orbits as ?_ — > oo. These solutions have 
a so-called isotropic singularity. 

CS a The fixed points CS a (a = 1,2,3) are hyperbolic with a three-dimensional unstable and a 
two-dimensional stable manifold; the former coincides with Bft , the latter is associated with 
the equations N^dNp/df- |cs Q = |(1 + 3u>) (for (3 ^ a). Therefore, each of the fixed points 
CS Q attracts an (equivalent) one-parameter set of type IX orbits in the limit f_ — > oo. 

TL Q Each fixed point of the line TL a has a three-dimensional center manifold and a two-dimen- 
sional stable manifold. The center manifold coincides with the vacuum subset Bft ft ; in 
Section 3.3 we have proved that the points of TL Q take the role of sources. We thus conclude 
that the fixed points on TL Q are center saddles. The two-dimensional stable manifold of 
each point of TL Q is contained in the LRS subset CRS a of B\x (which is the hyperplane 
given by the conditions Ep = S 7 and Np = iV 7 ); more specifically, the closure of the union 
of the unstable manifolds coincides with the closure of £IZS a . Therefore, for each fixed point 
on TL a there exists a one-parameter family of type IX orbits that converges to this point 
as f- — > oo; these orbits correspond to LRS solutions. (Conversely, generic LRS type IX 
solutions converge to TL a , see e.g. [1].) 

The solutions whose w-limit is one of the fixed points form a subfamily of measure zero of the 
(four-parameter) family of Bianchi type IX solutions. Following the nomenclature of [4] we thus 
refer to these solutions as non-generic solutions of Bianchi type IX. Alternatively, to capture the 
asymptotic behavior of these solution, we use the term past asymptotically self-similar solutions. 
(Since a fixed point in the Hubblc-normalizcd dynamical systems formulation corresponds to a 
self-similar solution, sec e.g. [1], solutions that converge to one of the above fixed points are 
asymptotically self-similar.) 

Apart from the invariant Bianchi contraction subsets there exists other invariant subsets of the full 
state space. The most important are the three equivalent LRS subsets CRS-y defined by E a = Ep 
and N a = Np for (af3j) e {(123), (231), (312)}. The past asymptotically self-similar solutions 
comprise the LRS Bianchi type IX solutions. As seen above, generic LRS solutions converge to 
TL Q towards the past (and each solution that converges to TL Q is LRS), but there exist exceptional 
LRS solutions that converge to F or CS a . The remaining orbits whose limit point is either F or CS Q 
correspond to past asymptotically self-similar solutions that are non-LRS. Clearly, every solution 
that converges to F or CS Q is a non- vacuum solution, since f2 ^ at F and CS Q . 

It is natural to ask how the non-generic orbits are embedded in the state space B\x- The LRS orbits 
form the three LRS subsets CRS a , which are the hyperplanes given by the conditions Ep = S 7 , 
Np = iV 7 , where (a/37) G {(123), (231), (312)}. The orbits whose w-limit set is the fixed point 
CS Q (for some a) form the set CS a in B\x; we call CS a the Collins-Stewart manifold. The local 
analysis of the fixed point CS a and the regularity of the dynamical system (8) imply that the 
Collins-Stewart manifold CS a is a two-dimensional surface; it can be viewed as a two-dimensional 
manifold with boundary embedded in Byx (where this boundary lies in _Bvn ). Analogously, the 
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orbits whose a- limit set is the fixed point F form the set 5" in B\x, which we call the isotropic 
singularity manifold, since solutions converging to F are those with an isotropic singularity. The 
isotropic singularity manifold 5" is a three-dimensional hypcrsurface; it can be viewed as a three- 
dimensional manifold with boundary. 

In the subsequent section we will state and prove the Bianchi type IX attractor theorem, which 
concerns the behavior of generic Bianchi type IX models (i.e., those that are not asymptotic self- 
similar, which provides an example of asymptotic self-similarity breaking; for other such examples, 
see [22]). 

5 A new proof of the Bianchi type IX attractor theorem 

Definition. Consider a solution of Bianchi type IX that is either vacuum or associated with a 
perfect fluid satisfying — | < w < 1. Such a solution is called generic if it is not past asymptotically 
self-similar. 

Remark. Accordingly, a solution is generic if its w-limit set is neither the point F, nor any of the 
points CS Q , nor a point on TL Q ; in other words, a generic solution corresponds to an orbit in Bix 
that is neither contained in the submanifold 5F, nor in CS a , nor in the hyperplane CRS a . Therefore, 
the set of generic Bianchi type IX states is an open set in B\x- 

Definition. The Mixmaster attractor Aix (alternatively referred to as the Bianchi type IX attrac- 
tor) is defined to be the subset of Byx given by union of the Bianchi type I and II vacuum subsets, 
i.e., Aix = -Bj ac ' U BJi°\ Accordingly, Aix consists of the three representations of the Bianchi 
type II vacuum subset and the Kasner circle (the Bianchi type I vacuum subset), i.e., 

A lx = K° U BT: U BT 2 U B V N ^ . (37) 

The main result concerning generic Bianchi type IX models is the Bianchi type IX attractor 
theorem, which is due to Ringstrom [4]; this result rests on earlier work that is reviewed and 
derived in [1] , and on [2] and [3] . In the following we state the Bianchi type IX attractor theorem 
in a version adapted to our purposes. 

Theorem 5.1 ([4]). A generic orbit 7 in B\x has an lj -limit set that is a subset of the Mixmaster 
attractor, i.e., C Aix = -Bf c - U B{f c - . 

Remark. Note that we have chosen the time direction towards the past singularity, see (6) and (13). 
When we use the standard future directed time, 'w-limit set' is replaced by 'a-limit set'. 

Remark. For an equivalent formulation of Theorem 5.1 let X(f-) = (£1, £2, £3, -^1,-^2, N3, H)(f-) 
be a generic solution of Bianchi type IX. Then 

\\X(f-)-Ant\\->0 (f_->oo), (38) 
where the distance \\X — A\x\\ is given as miny e _4 IX \\X — Y\\. 

Remark. Alternatively, the statement of the theorem can be expressed as follows: Along every 
generic orbit we have 

NiN 2 + N1N3 + N 2 N 3 -» and fi^O as f_ -> 00 . (39) 

In combination with the constraint (17b) we further obtain H — > 1 as r_ — > 00. Since, by 
definition, H = [l + ^(N 1 N 2 + NiA^ + N 2 N 3 ]~ 1/2 , sec (10) and (11), it follows that (N 1 N 2 + 
N1N3 + N2N3) when f_ — > 00. An immediate consequence, cf. (12), is that fi/Q — > 1 and 
dr-/d?- — > 1 as f— — > 00, so that r_ — > 00 as f— — > 00. Using these results, Theorem 5.1 instantly 
yields the original formulation of this theorem in [4]. 

In the following we give a new proof of Ringstrom's Bianchi type IX attractor theorem. The proof 
we present is subdivided into a number of lemmas, which culminate in Theorem 5.1. The additional 
remarks in this subsection are not needed directly for the proofs, but give further insights about 
the main ideas and the lines of argument. 
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Lemma 5.2. Every orbit in Bix possesses a non-empty co-limit set that is contained in the subset 
-Bviio of the boundary dB\x- 



Proof. Since the state space -Bix is relatively compact, every orbit has an w-limit point in Bix 
H U Byiio, cf. (20). The function 



A := H^N^Ns 
is positive and strictly monotonically decreasing on -Bix, since 

d 2 A 



dA 
dfl 



-3-1 A 
H 



dfi 



d 3 A 

s 2 =o dri 



£ 2 =0 



(40) 



2H- 1 ( s Sf + 3 S| + 3 S% ) A ; (41) 



in the perfect fluid case the first derivative is negative since q = 2£ 2 + i(l + 3w)Q > 0; in the 
vacuum case, i.e., O = 0, dA/df_ = is possible since q = when E 2 = 0, but then the third 
derivative is negative, since 3 5 2 + 3 5|+ 3 5| > when S 2 = because of the constraints. Application 
of the monotonicity principle to the function A yields that the w-limit set of every orbit is contained 
in dBix, where the set Hq is excluded because A = oo on Hq. □ 



The set -Bvn can be decomposed according to 

Eviio =BviioUBiiUBiU{Zi,Z 2 ,Z3}. (42) 

In the subsequent lemmas we exclude the possibility that generic orbits have cj-limits in Byilo or 
{Z Q }; for pedagogical reasons some of the lemmas refer to the (simpler) vacuum case and the fluid 
case separately. Furthermore, we prove that not only in the vacuum case but also in the fluid case 
only the vacuum subsets of B\ and B\\ come into question. Taken together these lemmas then 
directly lead to Theorem 5.1. 

Lemma 5.3 (Vacuum case). A generic orbit in B\x cannot have an ui -limit point on the subset 
-Bviio ■ 



Proof. {Vacuum case, i.e., 0, = 0. In this case 'generic' reduces to non-LRS, cf. Section 4.) In 
preparation for the proof we begin by considering a fixed point L a of the line TL a , see (25a); in 
particular, N a — 0, and Np = iV 7 > at L Q . The point L a has the following properties: L a is 
not hyperbolic, but on its center manifold, which is the set Byu Q1 there exists a one-parameter 
set of orbits converging to L Q as f_ — > — oo; in contrast, no orbit converges to L Q as f_ — ► oo; see 
the previous discussion of Bianchi type VIIo models. In this sense, L a resembles a (transversal) 
source for Bvn - Furthermore, from (15c) we obtain that N~ 1 dN a /df- = —6H at L a , hence 
L a possesses a one-dimensional stable subspace that lies in Bix- The center manifold reduction 
theorem applies [6]: There exists a neighborhood of L a such that the flow of the full nonlinear 
system is equivalent to the flow of the decoupled system 

—N a = -6HN a (43a) 
dr- 

-Jr-B VUo = F(Bvii ) , (43b) 

where i?vn denotes the collection of the variables of -Bvn ■ Since L Q is a (transversal) source for 
the second subsystem, L Q is a center saddle and there exists exactly one orbit whose w-limit is L a ; 
this orbit coincides with the unstable manifold of L Q (for which N a > 0, Np > 0, N-y > 0) and can 
straightforwardly be identified as a Bianchi type IX LRS orbit (because the vacuum type IX LRS 
subset is exactly solvable). 

Now consider a non-LRS orbit 7 in the vacuum subset of Bix and assume that the w-limit set 
^(7) contains a point on Byu , i.e., a point P such that N a = 0, Np > 0, N 7 > 0. We distinguish 
three possible cases: (i) P e TL Q ; (ii) P € B V n 0J p £ TL tt) P g QC a ; (iii) P £ QC a . 
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Consider case (i), i.e., assume P = L a e w(7). There are two possibilities: Either 01(7) = {L Q }, 
then 7 is the orbit that coincides with the stable manifold of L Q ; this is impossible since the stable 
manifold is an LRS orbit; or 0^(7) 3 {La}- The saddle structure of L Q allows us to draw the 
following conclusion: Since {L Q } C a; (7), it follows that also the stable manifold of L a must be a 
subset of (^(7). However, since the stable manifold lies in -Bix, this is a contradiction to Lemma 5.2. 
Hence, the existence of a point P G TL a in oj(j) is excluded for the non-LRS orbit 7. 

Case (ii) is analogous. Since P is contained in a; (7), so is the entire Bianchi type VIIo orbit through 
P and therefore the a- and w-limit sets of that orbit. The a-limit set is a fixed point h a on the 
line TL a , cf. the analysis of Bianchi type VIIo models; accordingly, L a D ^(7) which leads to a 
contradiction to Lemma 5.2 in analogy with case (i). 

Finally, consider case (iii). Since P G ^(7), the entire orbit through P (and the a- and w-Umit 
points of that orbit) must be contained in ^(7) as well; hence, QC a C a; (7). Because we have 
already excluded cases (i) and (ii), we know that £^(7) C BuU ( \J S QCs) ■ However, the heteroclinic 
orbits QCs, 5 = 1, 2, 3, are not connected, but 'isolated branches' of the set B\\ U ( [j s QCs) ■ Since 
such structures can never be part of a limit set, cf. the remark below, the assumption QC a C 0^(7) 
and therefore (iii) result in a contradiction. This finishes the proof of the lemma. □ 

Remark. Let us elaborate on case (iii) in the proof of the theorem and show explicitly why the 
isolated heteroclinic orbits QCs, 5 = 1,2, 3, are excluded from the possible w-limit set. In the proof 
of Lemma 5.3, case (iii), we assume that there exists an orbit 7 such that P G QC a is an element 
of w(7); hence there exists a diverging sequence of times, (<: n )neN, such that 7(^1) — > P (n — > 00). 
Let V be a sufficiently small neighborhood of Z Q , preferably generated by an open ball, so that 
P ^ V. For sufficiently large n, there exists times er„, <r„_i < a n < c„, such that j(a n ) G V. (This 
is a simple consequence of the continuous dependence of the flow on initial data; recall that QC a 
is a heteroclinic orbit connecting the fixed point Z a with the fixed point Q a .) Consequently, there 
exist times K n , S n -i < K n < S n , such the orbit 7 enters V at n n (i.e., 7(/t n ) G 9V, 7(re„ + e) G V 
for sufficiently small e > 0). For all n, 7(/t„) is contained in the complement of a sufficiently small 
neighborhood 11 of the point dV n QC a (because the flow of the dynamical system points out of 
V in dV fl It). By going over to a subsequence 7(k„), this implies that 7 has a w-limit point on 
dV\U, i.e., an w-limit point in -Bvn (or Ho) that is not contained on QCs, 6 = 1,2,3. This is a 
contradiction to the assumption oj(j) C B\\ U (|J 5 Q£i). 

Remark. The concept of 'isolated branches' — used in the proof of Lemma 5.3 and discussed in the 
previous remark — is a very useful picture to have in mind also for the general situation. The state 
space Bix can be depicted roughly as the space between three branches (the three -Bvn subsets 
^NtNo ' Bn 2 n 3 -> ^n 3 n ± ) sticking out from a common basis ('trunk') which is the set B\\; the end 
points of the branches are the points Zi, Z2, Z3; see Figure 3. The main observation is that the flow 
on the branches is (more or less) unidirectional, i.e., directed away from the points Z a (although 
there are the lines TL a of fixed points on JBvn which make the situation much more subtle). 
A generic orbit 7 cannot have w-limit points in the interior of the state space, see Lemma 5.2, 
however, the picture strongly suggests that w-limit points on the branches arc excluded as well: 
This is simply because an orbit 7 cannot 'climb up' to a point on -Bvti along the -Bvii branches 
themselves (which is due to the continuous dependence of the flow on initial data). Lemma 5.4 is 
intimately connected with this idea. 

Remark. Lemma 5.3 implies that the w-limit set of a non-generic orbit 7 (in the vacuum subset of 
Bix) is a subset of Bn U {Z Q }. It follows trivially that {Z a } C 0^(7) is impossible, which is simply 
because these points are isolated from the set B\\ (while limit sets are necessarily connected). (By 
studying the LRS subset it is easy to show that {Z Q } C 01(7) is excluded for LRS orbits.) We will 
give an independent proof of this statement (and its extension to the fluid case) in Lemma 5.4. 

Lemma 5.4. Let 7 be an orbit in B\x- Then 0^(7) D {Zi} (or {Z2}, {Z3}) is impossible. 

Proof. Since the statement of the lemma is trivial for non-generic orbits (past asymptotically self- 
similar orbits, cf. Section 4), we restrict ourselves to generic orbits 7. Assume that there exists a 
(generic) orbit 7 such that the w-limit set of 7 contains Z 1; i.e., 0^(7) 3 {Zi}. Accordingly, there 
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Figure 3: A simplistic picture of the state space -Bix- Bvn consists of three 'branches' (Bft ft , 
^jViAt 3 i Bft 2 ft 3 ) that intersect in a common subset, B\\. The flow of the -Bvn subsets is directed 
away from the points Z Q towards Bu. 



exists a diverging sequence of times (<r n )neN such that 7(^1) — > Zi as n — > oo. There exists a 
neighborhood V of Zi such that 7 intersects the complement of V infinitely many times (otherwise 
7 would converge to Zi). Therefore we can construct a sequence of times (cr n ) n eN, o~ n < <Zn Vn, such 
that 7(cr„) € 9V and 7|( CT „.?„] G f° r au n - (Clearly, |c n — cr„| diverges as 71 — > 00, because 7(? n ) 
converges to the fixed point Zi.) The sequence (7(cn)) n gN possesses a converging subsequence, i.e., 
there exists a point P such that 7(cr n ) — > P as n — > 00 (where the index n now runs over the index 
set of the considered subsequence). By construction, P G dV; furthermore, by definition, P is an 
w-limit point of the orbit 7; Lemma 5.2 implies that P G -Bvn (since B\ and Bu are disjoint from 
V). As discussed in the subsection "Bianchi type -Bvn : A new analysis," the w-limit of the orbit 
through the point P, which we call 7p, is a fixed point on TLi or Bu (gcnerically, w(tp) is a fixed 
point on the Kasncr circle); in particular, w(tp) does not contain Zi. Assume that w(7p) G -Bu; 
then, by continuous dependence on initial data, for sufficiently large n, 7|( (T , l! .) shadows the orbit 
7p and reaches Bu; but this is a contradiction to the fact that 7| („•„,.) C V. Assume that oj(jp) 
is one of the fixed points of TLi ; if this fixed point is not contained in V, we immediately obtain 
a contradiction; otherwise, we proceed in close analogy to the proof of Lemma 5.3; we exploit the 
(center) saddle property of w(7p) and we obtain that 7|(o- n .), for sufficiently large n, shadows first 
7p, then follows some orbit emanating from u>(yp) and eventually approaches a fixed point on the 
Kasner circle; this is again the desired contradiction. □ 

Lemma 5.5 (Fluid case). A generic orbit in Bix cannot have an u -limit point on the subset 
Bvu ■ 

Proof. [Fluid case, i.e., Cl > 0.) Consider a generic orbit 7 in _Bjx and assume that the w-limit 
set a; (7) contains a point P on Byu . We distinguish three possible cases: (i) P is an element of 
the fluid subset of -Bvn , i.e., fi|p > 0; (ii) P G QC a ; (iii) P G TL Q ; (iv) P is an element of the 
vacuum subset of Bvn 1 but P TL Q , P QC a . 

Consider case (i). Since P G ^(7), the orbit 7p through P and its a-limit cc(7p) must also be 
contained in Our analysis of the -Bvn subset shows that a(7p) = {Z a }. Consequently, 

{P,Z a } C 0^(7); but this is a contradiction to Lemma 5.4. Case (ii) is completely analogous. 

Consider case (iii), i.e., 0^(7) contains a fixed point L Q of the line TL Q . The scenario 0^(7) = {Lq,} 
is impossible, since 7 is a generic orbit and the stable manifold of L Q is a subset of the LRS subset. 
In analogy to the proof of Lemma 5.3 we can exploit the (center) saddle structure of L Q , which 
is reflected in (43) and the additional 'fluid' equation ti~ 1 dti/df_ — —3(1 — w)H. Hence, since 
{L Q } C uj(j), it follows that 0^(7) contains a point Q on the stable manifold of L a ; since Q G" Bix 
by Lemma 5.2, we have Q G -Byiio ^ CRS a , i.e., Q lies on the Bianchi type VILj LRS orbit that 
converges to L a . Because Cl > for this orbit (and thus for Q), this brings us back to case (i); a 
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contradiction ensues. 

Consider case (iv). Since the a-limit of the orbit through P is a fixed point L Q on TL a , case (iv) 
can be reduced to case (hi). □ 

Lemma 5.6. A generic orbit in Bix cannot have an ui -limit point on the fluid subset of B\\. 



Proof. The lemma is obviously true in the vacuum case. In the fluid case, assume that there exists 
a generic orbit 7 such that ^(7) contains a point PeBiU JE?n with f2|p > 0. First, assume that 
P is an element of Bu, e.g., in the B$ a subset. The orbit through P has the fixed point CS Q as 
its a-limit set, see Section 3; hence CS a £ ^(7). Using the saddle structure of CS a , see Section 4, 
we conclude that 07(7) also contains a point Q of the stable manifold of CS Q . However, since this 
stable manifold is a subset of -Bvn U -Bix, we have Q £ -Bvn U -Bix and thus a contradiction to 
Lemma 5.2 or Lemma 5.5. Second, assume P € By. Using the same line of arguments (where F 
takes the role of CS a ) we obtain that Q £ Bu U -Bvn U -Bix and thus a contradiction to what has 
already been proved. □ 



Lemmas 5.3 and 5.5 imply that the cj-limit set of a generic orbit in Bix must be either contained 
in the vacuum subset of Bu (where we recall that Bu = Bu UBi) or it is a fixed point of the set 
{Zi, Z2, Z3}. It remains to prove that the latter scenario is impossible. 

Lemma 5.7. There does not exist any orbit in Bix that converges to Zi (L 2 , L3) as f_ — » 00. 



Proof. The non-trivial case (which is at the same time the case that is relevant for our purposes) is 
the non-LRS case. We perform a proof by contradiction. The main idea is to consider a non-LRS 
orbit that is assumed to converge to the fixed point Zi; the convergence to Zi then implies that N± 
decays rapidly so that the orbit is forced to shadow a -Bvn orbit as Zi is approached. However, 
there docs not exist any orbit on -Bvn that converges to Zi as f_ — ■> 00; a contradiction must 
ensue. 

This idea is formalized by using a non-negative function £1 on Bix which is zero at Zi and whose 
restriction to JBvn is monotonically increasing along the flow of Byu . Since the orbit 7 in JBix 
that is assumed to converge to Zi must shadow a type VIIo orbit, we expect the function £1 to 
increase along 7 as well, or at least to decrease at a rate small enough so that the integral still 
exists. This leads directly to a contradiction, because by assumption £1 must go to zero as Zi is 
approached. In essence this is the main idea employed by Ringstrom [4] in a similar context, who, 
however, used ii/-normalized variables, and a function introduced by Wainwright and Hsu [23]; to 
facilitate comparison we choose basically the same function. Consider 

c, - (£ *- £ ' ) ; + J "> 2 , (44,) 

tit = ~N^N~ 3 t 2( ^ 2 " ^ 3)2( ^ 1 ~ 2ll) + 2{t ' 2 ~ ^ )( ^ 2 ~ • (44b) 

The second term in the brackets can be simply estimated by using 

2(E 2 - E 3 )(iV 2 - N 3 ) < (£ 2 - £ 3 ) 2 + (N 2 - N 3 ) 2 = bN 2 N 3 . (45) 
We combine the constraints (17a) and (17b) to find 

Y?-H 2 + ± [Nl + (N 2 - N 3 ) 2 ] + n = iJVjfJVa + N 3 ) . 
Since S 2 = ill 2 + ±(E 2 - S 3 ) 2 we obtain 

i(S x - 2#)(£ 1 + 2H) + i [iV 2 + (E a - E 3 ) 2 + (N 2 - N 3 ) 2 ] + Q = ±Ni(N 2 + N 3 ) . 
If Ex — 2H > 0, which is the 'worst case scenario' for our considerations, then 

1--- 4 1 _ - TVi Ni 

^l-2H< -Ni(N 2 + iV 3 ) - i + 2J? < -(N 2 + N 3 )^<\^ (46) 
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for some sufficiently large positive constant A (since N a are bounded); if Si — 2H < 0, this 
inequality holds trivially. 

Assume that there exists a non-LRS orbit 7 in Bix that converges to Z\ as f_ — > 00. Inserting 
the inequalities (45) and (46) into (44b) yields 



(47) 



where the latter inequality is true for sufficiently large r_, since H — > along 7. At a reference 
time f_ = f , which we choose to be sufficiently large (in order for (47) to hold), the function Ci 
takes a value £1 > 0. Integration of the differential inequality (47) yields 



(48) 



along the orbit 7. 

To estimate the r.h. side in equation (48) we exploit the convergence of 7 to Z\. This convergence 
entails that Ni — -> and N2 — ► \/6, A^3 — > v 7 ^; the constraints then automatically imply H — > 0, 
£ 2 — > 0, and f2 — > 0. Let <5 = 2V6 — N2 — Ns- Convergence to Zi implies that <5 — > 0. Using the 
constraint (17b) and the fact that Ni = 0(H 3 ) as H — > (which follows from the monotonicity of 
A) we obtain 



H 



7E S 



o(5) 



H 2 



7E S 



0(6) 



for 8^0. 



(49) 



The constraint (17a) gives 



n = -M 



o(<5) 



H 2 



= 1 + o(l) for S 



0. 



(50) 



For the quantity q = 2(S 2 + f2) — |(1 — w)£l we therefore derive the estimate q/H 2 = 2 + o(l) — 
|(1 - w)Ct/H 2 > i(l + 3w) + o(l), while g/-ff 2 = 2 + o(l) in the vacuum case. We conclude that 
there exists a constant a > such that 

- ;; _ <-> 

for sufficiently small values of 5, i.e., for sufficiently large values of t_ along 7. Using this inequality 
in the integrated version of equation (41), i.e., 



A 



N1N2N3 
H 3 



oc exp 



-3 



tdf. 
H 



(52) 



results in the estimate 



-p<Cexp 



Hdf'_ 



(53) 



where C is some positive constant and where r_ is assumed to be sufficiently large. 

Therefore, N\ goes to zero at a fast rate which ensure finiteness of the integral in (48) as f_ — > 00: 



iVi 



Ni 



^4df_ = I ^Hdf_ < GaT 1 



1 — exp —a 



Hd?_ 



< Ca 



(54) 



We thus conclude from (48) that £1 remains bounded away from zero as f_ — * 00; but this 
contradicts the assumption that the orbit 7 converges to the fixed point Zi, since then £1 — * 
along 7. □ 

Remark. The function £1 can be used to give an independent proof of the statement Zi <f_ Lu("f) V7 
(which is closely related to Lemma 5.4); we briefly sketch this proof. Since the case ^(7) = {Z\} is 
treated in Lemma 5.7, consider an orbit 7 such that 01(7) C {Zi}. If 7(<r n ) ~~ ► Zi as n — ► 00, then 
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Ci(^n) — > as 7i — > oo. Relations (47) and (51) hold whenever an orbit is in a sufficiently small 
neighborhood of Zi. Therefore, in analogy to the considerations in the proof of Lemma 5.7, it is 
impossible to achieve Ci(<Tn) — * 0, if 7 has an w-limit point that does not lie on the LRS subset. 
(In there were such a point it would be impossible for £1 to decrease sufficiently much between 
that point and 7($»)j n sufficiently large.) A priori it is possible that 7 has an w-limit point on 
the LRS subset (note that (1 = on the LRS subset). However, a study of the LRS dynamics 
shows that orbits emanate from Zi, but there do not exist orbits that converge to Zi as f_ — > 00. 
Using the same reasoning as above (e.g., in the remark following Lemma 5.3 or in the proof of 
Lemma 5.4) we therefore exclude the possibility that 7 has an w-limit point on the LRS subset as 
well. Consequently, ^(7) cannot have an w-limit point except Z^ but this is a contradiction to the 
assumption. 

The collection of the lemmas finally yields Theorem 5.1 and thus completes our argument. 

Theorem 5.1. A generic orbit 7 in B\x has an co-limit set that is a subset of the Mixmaster 
attractor, i.e., a; (7) C Aix = Bf- C - U £?j v f c ' • 

Proof. Lemma 5.2 implies that the cj-limit set of a generic orbit in Bix is a subset of -Bvn 1 which 
is given by (42). Lemmas 5.3 and 5.5 exclude -Bvn , Lemma 5.6 excludes the fluid subsets of 
B\ and B\\. Lemma 5.7 (in combination with Lemma 5.4) excludes {Zi,Z2,Z3}. This leaves the 
vacuum subset of B\ U Bu as the only possible superset of the w-limit set of a generic orbit. □ 

Remark. To complete the statement of Theorem 5.1 it is important to note that both ^(7) C B j ac ' 
and a; (7) C BJf c ' is impossible. The latter is obvious, since BJf c ^ consists of a collection of 
heteroclinic orbits (transitions) with end points on £?j ac - (= K°), see Section 3. The proof that 
^(7) Q B\ &c - (= K°) is impossible is contained in the local analysis of Section 4. (For an alternative 
proof see [3].) 

Remark. It is important to emphasize that the Theorem states that (^(7) C Aix- Whether a; (7) 
actually coincides with Aix (at least generically) or whether it is a proper subset of Aix is open. 
This question and related issues are discussed in detail in [10]. 



6 Consequences 

The Bianchi type IX attractor theorem in conjunction with our understanding of the flow on the 
attractor subset implies a number of consequences that we formulate as corollaries. Some of these 
correspond to results presented in [4]; our approach, however, is rather different. 

On the Mixmaster attractor Aix, the dynamical system (15) generates an intricate network of 
structures that are invariant under the flow: heteroclinic cycles and finite and infinite heteroclinic 
chains. (These heteroclinic structures arise by concatenating the vacuum type II orbits on , 
Bfj 2 , and Bfj 3 , see Figure 2.) The 'simplest' structure (i.e., the structure that contains the smallest 
number of fixed points) is a heteroclinic cycle with three fixed points on K°, sec Figure 4. We 
refer to [10] and references therein for a comprehensive discussion. 

Lemma 6.1. If P £ Aix is an lu -limit point of a type IX orbit, then the entire heteroclinic 
cycle/chain through P must be contained in the u -limit set. 

Proof. The lemma follows from basic facts of the theory of dynamical systems [5]. □ 

Lemma 6.1 entails that the heteroclinic cycles and chains are potential limit set candidates for 
generic type IX orbits. 

Corollary 6.2. The asymptotic behavior of a generic Bianchi type IX solution is oscillatory with 
oscillations between at least three fixed points on K°. 
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Figure 4: The 'simplest' structure on the Mixmaster attractor: a heteroclinic cycle with three 
fixed points. In general, concatenating of the vacuum type II orbits of Figure 2 yields infinite 
heteroclinic chains. 

Proof. The simplest structure on Aix that is a potential cj-limit set for a generic type IX orbit is a 
heteroclinic cycle with three fixed points on K°, see Figure 4. A type IX orbit converging to such 
a heteroclinic cycle (if such an orbit exists) exhibits oscillations between three Kasner points. □ 

Remark. The oscillatory behavior implies that generic asymptotic type IX dynamics constitute an 
example of asymptotic self-similarity breaking [22]. 

Corollary 6.3. If one of the Taub points {Ti, T2, T3} is an u> -limit point of a type IX orbit, then 
the uj -limit set contains Kasner fixed points arbitrarily close to the Taub points. 

Proof. Assume the contrary, i.e., suppose that there exists a type IX orbit 7 such that T a (for 
some a) is an element of ^(7), while at the same time there exists a neighborhood U of T Q such 
that w(7)nWnK° = {T a }. However, uj(-f)C\U 3 {T Q }, since £^(7) is connected and strictly larger 
than T a . (There do no exist type IX orbits that converge to a Taub point, see Section 4.) Taking 
into account the structure of orbits on Aix we conclude that a; (7) = ^(IW, where ? denotes 
the closure of the type II orbit Q a — > T Q . However, since w-limit sets cannot contain 'isolated 
branches' of this type, we obtain a contradiction. □ 

Remark. Corollary 6.3 implies that the cj-limit set contains an infinite set of Kasner fixed points 
in a neighborhood of the Taub point(s), but this set is not necessarily a continuum of fixed points. 

Lemma 6.4. Let W be a neighborhood of the Mixmaster attractor Aix- Then there exists a smaller 
neighborhood V of Aix, V C W, such that each solution with initial data inV (at f- = 0) remains 
in W for all r_ > 0. 

Proof. Henceforth, we denote by X(x;f-) the type IX solution generated by initial data x G 
Six, i.e., X(x,f-) is a type IX solution and X(x,0) = x. We employ the two functions and 
N1N2 + N1N3 + N2N3, cf. (39), as a measure of the distance from the Mixmaster attractor. 
The choice of a small neighborhood of the Mixmaster attractor then corresponds to both and 
N1N2 + NxNz + N 2 N 3 being bounded by a small constant e. Choose a small neighborhood W of 
Aix (corresponding to a choice eyy) and consider a neighborhood V of Aix with V C W, ey -C eyy- 
A solution X(x,f-) with initial data x £ V remains in W at least for < f_ < TL. The constant 
T_ depends on ey and, a priori, on the initial data; we intend to show that T_ = 00 irrespective 
of the choice of initial data if ey is sufficiently small. 

First, consider Eq. (18) for O. At the Kasner circle we obtain rt~ 1 dti/d?^ | K o = —3(1 — w); 
hence there exists a neighborhood U\ of K° such that Vt is (exponentially) decreasing as long as 
the solution stays U\. In general, we obtain VL~ 1 dVL/d?^ < (1 + 2>w) in VV, i.e., a limit on the 
possible increase. Decompose W into Wi = W n U\ and W11 = W\U\. Using the flow on _4dc 
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(and elementary results from the theory of dynamical systems) we see that the solution X(x, f_) 
oscillates between Wi and Wn, where the sojourn times in Wi are large compared to the sojourn 
times in VVn. (The ratio of the respective sojourn times diverges with ey — > 0.) Using the estimates 
for dti/dr- it follows that f2 exhibits a (rapid) overall decrease for f_ £ [0,T-) (i.e., as long as 
X(x,f-) £ W). Consequently, at r_ = TL, the quantity O is smaller than initially, £l\f_ < ey. 

Second, consider_7V := N-lNs + N1N3 + N2N3. For (a/37) £ {(123), (231), (312)} wc obtain 
(NaNp)' 1 d(N a Np)/d?_ =_-2(qH - £ 7 +_F) from (15c). Let W T be a (small) neighborhood 
of the Taub points. Then N a N/3 and thus N is (exponentially) decreasing in Wi\Wr- In Wn we 
obtain N~ 1 dN /df_ < const, i.e., a limit on the possible increase. To obtain information about TV 
in Wt we use the center manifold reduction theorem: The set Wt contains a piece of the Taub line 
TL 7 (up to some value of N a = Np). The analysis of -Byiio in Section 3.3 implies that each fixed 
point L 7 on TL 7 is the (non-hyperbolic) source for a one-parameter family of orbits in -B V n ■ The 
family of orbits emerging from L 7 forms a two-dimensional invariant surface <S[L 7 ] in Byn H Wt- 
The proof of Lemma 3.1 entails that there exists a universal constant ct such that 

sup 5[ L 7 ] (N a Np) < c T inf s w (N a Np) (55) 

for each surface <S[L 7 ] (in fact, ct is only marginally larger than 1). We thus obtain, in D Wt, 

a foliation of -Byji m ^o two-dimensional invariant surfaces across which the relative variation of 
(N a Np) is small. Since each point L 7 is a center saddle with a two-dimensional stable sub- 
space/manifold, the foliation of -Byn C carries over, by the center manifold theorem, to a 
foliation of Wt into four-dimensional invariant hypersurfaces 7i[L 7 ]. (Each hypersurface is associ- 
ated with the direct sum of the stable subspace and the tangent space of <S[L 7 ] at L 7 .) Furthermore, 
for sufficiently small Wt, Eq. (55) carries over from <S[L 7 ] to 7i[L 7 ]. Therefore, since every solution 
in Wt is contained in one of the hypersurfaces 7i[L 7 ], we obtain (N a Np)(f-) < CT(N a Np)(a-) 
for all t_ such that the solution is contained in Wt, where cr_ is the time the solution enters Wt- 
(Note in particular that this result is independent of the actual time the solution spends in the 
neighborhood Wt-) Summing up: For the solution X(x, r_) we find a decrease of N in Wi\Wt 
(comparatively long sojourn times), an increase in Wn (comparatively short sojourn times) and 
a bound on the increase in Wt (irrespective of the sojourn times). From these facts we infer an 
overall decrease of N for f_ £ [0, TL) (i.e., as long as X(x,f_) £ W). Consequently, at f_ = T_, 
the quantity N is smaller than initially, N\f_ < £y- 

Since both < ey < eyy and N\f_ < e v < ew, the solution X(x, f_) does not actually leave W 
at f_ but is still in V (c W) for some time beyond TL. It is immediate that this leads to T_ = 00, 
and hence the lemma is established. □ 

Remark. Note that not all solutions in V or W are generic. In these neighborhoods of Aix, there 
also exist LRS solutions, which converge to TL Q , a = 1, 2, 3, instead of to Aix- 

Corollary 6.5. Convergence to the Mixmaster attractor is uniform on compact sets of generic 
initial data: Let X be a compact set in Bix that does not intersect any of the manifolds 5", CS a , 
CRS a , so that each initial data x £ X generates a generic type IX solution. Let X(x;f-) denote 
the type IX solution with X(x, 0) = x. Then 

||X(^f_)-^ IX |H0 (f_->oo) (56) 

uniformly in x £ X. 



Proof. Let W be a neighborhood of the Mixmaster attractor Aix- We need to show that there 
exists T- such that X (&;?-) £ W for all r_ > TL, for all x £ X. Let V be as in Lemma 6.4. 
Theorem 5.1 implies that for each x £ X there exists &-{x) such that X(x; f_) £ V. Since X is a 
compact set of generic initial data, XL := sup £G3C P-(x) exists. (If this supremum did not exist, we 
could consider a sequence of initial data along which <r_ diverges to construct generic initial data 
violating Theorem 5.1.) Applying Lemma 6.4 we conclude that X(£; f_) £ W for all f_ > TL, for 
all x £ X. □ 
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Corollary 6.6. For generic solutions of Bianchi type IX the Weyl curvature scalar C abc dC 
(and therefore also the Kretschmann scalar) becomes unbounded towards the past. 

Proof. At a fixed point on K° (which represents a Kasner solution) the Hubble-normalized Weyl 
curvature scalar C abcd C abcd / '(48iJ 4 ) is given by 27(2 - SiS 2 S 3 )- Therefore, C abcd C abcd / '(48iJ 4 ) € 
[0,4] on the Kasner circle K°, where C abcd C abcd / (48ff 4 ) = holds only at the Taub points 
{Ti,T 2 ,T 3 }. Since the w-limit set of a generic type IX solution must necessarily contain a fixed 
point on K° different from the Taub points by Corollary 6.3, we conclude that C a bcdC abcd becomes 
unbounded towards the past; we simply use that H — > oo as f_ — > oo (t — > 0). □ 

Corollary 6.7. Taking into account both the expanding and contracting phases of Bianchi type IX 
solutions, generic Bianchi type IX initial data generate an inextendible maximally globally hy- 
perbolic development associated with past and future singularities where the curvature becomes 
unbounded. 

Remark. This is a direct consequence of the previous corollary. It follows straightforwardly that the 
analogous statement holds for the asymptotically self-similar solutions as well, the only exceptions 
being the type IX vacuum LRS solutions. 

7 Discussion 

In this paper we give new and comparatively short proofs of the main rigorous results on Bianchi 
type IX asymptotic dynamics: Ringstrom's Bianchi type IX attractor theorem [4], Theorem 5.1, 
and its consequences. To find more succinct arguments is not our primary motivation to re- 
investigate the problem. By emphasizing the importance of the Lie contraction hierarchy our 
proof demonstrates that Bianchi type IX is special in comparison with the other oscillatory Bianchi 
types: type VIII and type VI_i/ 9 . Let us elaborate. 

Among the class A Bianchi models Bianchi type IX is characterized by the condition that the three 
structure constants possess the same sign, see Table 1; in terms of the dynamical systems variables 
of Section 2 positivity of the structure constants is expressed by the conditions Ni > 0, N2 > 0, 
N3 > 0. The pivotal feature of type IX dynamics is the following fact: The set of asymptotic 
states that are accessible to Bianchi type IX models in the past asymptotic limit is represented 
by the 'Bianchi type IX Lie contraction hierarchy' of Figure 1. Since this hierarchy is obtained 
by successively setting the type IX structure constants to zero (Lie contractions), the first level is 
taken by the three equivalent representations of the Bianchi type VIIo state space, the second level 
consists of the three representation of Bianchi type II, and the third level coincides with Bianchi 
type I, see Figure 1. Accordingly, admissible past asymptotic states of type IX models are of 
Bianchi type I, where JVi = 0, N 2 = 0, N 3 = 0, of Bianchi type II, where N a > 0, Np = A 7 = 0, 
or of Bianchi type VIIo, where N a > 0, Np > 0, A 7 = 0; cf. Figure 1. (As usual, (aflj) runs 
over the set {(123), (231), (312)}.) Models of Bianchi type VIo are not among the admissible past 
asymptotic states, since type VIo does not appear in the Bianchi type IX Lie contraction hierarchy; 
this is because one of the structure constants is necessarily negative in type VIo- (Using dynamical 
systems terminology to summarize: The w-limit set of every type IX orbit lies on the boundary of 
the type IX state space, which is the union of the type I, type II, and type VIIo state spaces; in 
contrast, the type VIo state space is not part of the boundary.) 

For simplicity and clarity, let us restrict our discussion to the vacuum case. Vacuum Bianchi type I 
solutions (Kasner solutions) are represented by fixed points on the Kasner circle, K°, while vacuum 
Bianchi type II solutions are represented by heteroclinic orbits that connect one fixed point on K° 
with another; see Figure 2. The Mixmaster attractor (which is simply the union of the Kasner circle 
and the three equivalent representations of the type II vacuum subset) is covered by an intricate 
network of these heteroclinic orbits, which can be concatenated to form heteroclinic cycles and 
chains. While Bianchi type I and type II thus fit together seamlessly to form the fabric of type IX 
asymptotics, Bianchi type VIIo is the odd one out. Solutions of type VIIo do not connect one fixed 
point on the Kasner circle with another, but connect a fixed point on K° with a fixed point on the 
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Taub line, for which N a = Np > 0. As a consequence, Bianchi type VIIo orbits are incompatible 
with the network of heteroclinic cycles/chains; in particular, Bianchi type VIIo orbits cannot be 
concatenated with Bianchi type II orbits to form heteroclinic chains. This incompatibility is crucial: 
It can be regarded as the underlying reason that excludes Bianchi type VIIo from the set of possible 
past asymptotic states. Indeed, the new proof of Theorem 5.1 given in this paper reflects this idea 
accurately. 

The exclusion of Bianchi type VIIo from the set of admissible past asymptotic states is the cor- 
nerstone of the analysis of type IX asymptotics. Once established, the exclusion of type VIIo 
automatically leaves the Mixmaster attractor (or a subset thereof) as the past attractor and thus 
yields Theorem 5.1. So what about Bianchi type VIII then? There is but one difference be- 
tween type VIII and type IX that is relevant in the context of past asymptotic dynamics: The 
Lie contraction hierarchies of Bianchi type VIII and Bianchi type IX differ on the first level, com- 
pare Figure 1 with Figure 5. The boundary of the type VIII state space encompasses a Bianchi 
type VIIo state space and two equivalent representations of the Bianchi type VIo state space in- 
stead of three type VIIo representations; hence the Bianchi type VIII Lie contraction hierarchy 
contains a representation of each class A model, whereas Bianchi type IX does not. (This is merely 
one aspect of Bianchi type VIII being more general than type IX. Another aspect is the violation 
of the permutation symmetry in type VIII.) 

D 

-Byill 5 



-Bvi /-Bvii 4 



Bu 3 



Bj 2 

Figure 5: This subset contraction diagram for Bianchi type VIII is to be compared with the one 
of type IX, see Figure 1. The subscript of B* denotes the non-zero variables; variables with the 
superscript + are are positive, variables with the superscript _ are negative; e.g. B^-^+ denotes 
the type VIo subset with N\ < 0, V2 > and V3 = 0. Note that B^- is a completely equivalent 
representation of B^+. 

Bianchi type VIo and Bianchi type VIIo possess the same scale-automorphism group which gen- 
erates corresponding monotone functions [1]. However, the state spaces on which these monotone 
functions act are rather different: While the type VIIo state space is unbounded and contains a 
line of fixed points, the type VIo state space is bounded and does not contain this line of fixed 
points. The reason for these differences between type VIo and VIIo is that type VIo does not admit 
an LRS subset, while VIIo docs. In type VIIo the vacuum LRS subset contains a one-parameter 
set of representations of Minkowski spacetime which are different, but equivalent, from the Taub 
representation in Bianchi type I. In the dynamical systems picture this corresponds to a line of 
fixed points, TL Q , which is necessarily absent in type VIo- These features indicate that type VIo 
dynamics is simpler than that of type VIIo, which is indeed the case, see, e.g., [9, 1, 21]. Paradox- 
ically, the very simplicity of the type VIo flow is the reason why the analysis of Bianchi type VIII 
asymptotics is more intricate and complex than that of type IX. Let us elaborate. In contrast to 
Bianchi type VIIo solutions, orbits of Bianchi type VIq are compatible with the network of hctc- 
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roclinic chains on the Mixmaster attractor. This is a crucial fact, because it implies that Bianchi 
type VIo orbits can be concatenated with heteroclinic chains. Therefore, Bianchi type VIo states 
cannot be ruled out a priori as possible asymptotic states of Bianchi type VIII solutions. It is 
clear that this causes the (analog of the) proofs of Theorem 5.1 to fail in type VIII. However, the 
failure of the methods of proof does not imply the failure of the statement: It is expected that (a 
generic version of) Theorem 5.1 holds and that the Mixmaster attractor is the past attractor for 
type VIII models. However, it is conceivable that arguments of a completely different kind, like 
stochastic arguments, are necessary to prove that Bianchi type VIo is excluded (generically) from 
being involved in the asymptotic dynamics of solutions. A statement like this would then be the 
core of a (generic) version of Theorem 5.1 for Bianchi type VIII. 

A closely related issue concerns the validity of Lemma 6.4 and the related Corollary 6.5. Lemma 6.4 
fails in Bianchi type VIII. There exists (generic) initial data that is arbitrarily close to the Mix- 
master attractor but transported beyond a given neighborhood of the attractor by the flow of 
the dynamical system. The reason for this difference between type VIII and IX is the by now 
familiar one: The contraction hierarchy of type VIII contains type VIo while that of type IX does 
not. However, lack of uniform convergence does not imply that there is no convergence at all in 
type VIII: The Mixmaster attractor might still be the past attractor. 

Our arguments support the thesis that Bianchi type IX is special. The existence of discrete symme- 
tries associated with axes permutations, which follow from the positivity of the structure constants, 
leads to a simplification of the problem and makes the treatment of type IX dynamics relatively 
straightforward. Type VIII on the other hand is less symmetric and is thus harder to grasp. 
The same is true for the other remaining oscillatory Bianchi models — the Bianchi type VI_!/g 
models. These types might therefore be more relevant for our understanding of generic spacclikc 
singularities. 
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